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Abstract. The term active nematics designates systems in which apolar elongated 
particles spend energy to move randomly along their axis and interact by inelastic 
collisions in the presence of noise. Starting from a simple Vicsek-style model for 
active nematics, we derive a mesoscopic theory, complete with effective multiplicative 
noise terms, using a combination of kinetic theory and Ito calculus approaches. 
The stochastic partial differential equations thus obtained are shown to recover 
the key terms argued in EPL 62 (2003) 196 to be at the origin of anomalous 
number fluctuations and long-range correlations. Their deterministic part is studied 
analytically, and is shown to give rise to the long-wavelength instability at onset of 
nematic order (see arXiv:1011.5408). The corresponding nonlinear density-segregated 
band solution is given in a closed form. 
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1. Introduction 

The study of collective properties of systems of interacting active particles [H [21 E] is 
currently attracting a great deal of interest. In active matter, particles extract energy 
from their surrounding and dissipate it to propel themselves in some coherent way 
in a viscous fluid and/or over a dissipative substrate. In this last case, or whenever 
hydrodynamic effects can be neglected, physicists speak of "dry active matter" [3]. 
Systems as diverse as animal flocks [H El [6] , human crowds [71 [8] , subcellular proteins 
[9], bacterial colonies [10], and driven granular matter [Til 1121 [13] have been described 
in this framework. 

In the context of dry active matter, there is now some consensus in the physics 
community that minimal models such as the celebrated Vicsek model [IH [IS] play a 
crucial role, since they stand as simple representatives of universality classes which 
have started to emerge from a combination of numerical and theoretical results: for 
instance, many different microscopic (particle) models have been shown to exhibit the 
same collective properties as the Vicsek model, and the continuous equation proposed 
by Toner and Tu [TB] is widely beheved to account for its collective properties. Such 
hydrodynamic theories formulated at the mesoscopic level (stochastic PDEs) are the 
natural framework to characterize and deflne universality classes. 

In early approaches these mesoscopic theories have been built on the principle 
of including all that is not explicitly forbidden, retaining all leading terms (in a 
gradient expansion sense) allowed by symmetries and conservation laws [IHITT]. This 
grants access to the general structure of these equations and has been successful in 
describing relevant features of active matter systems such as their anomalously large 
number density fluctuations [161 [121 IlHl [IS]. Despite the attractions of a gradient 
expansion, it typically contains many transport coefficients of unknown dependence on 
microscopic control parameters and hydrodynamic fields such as local density. Moreover, 
the dependence of the noise terms on the dynamical fields in such equations remains 
arbitrary, and frequently neglected, whereas it could have profound consequences 
for important phenomena such as spontaneous segregation, clustering and interface 
dynamics. 

Ideally, thus, one would be able to derive well-behaved mesoscopic theories using 
a systematic procedure starting from a given microscopic model. Kinetic-theory-like 
approaches [13 [201 Ell [221 [23] go one step towards this goal, by allowing one to 
compute hydrodynamic transport coefficients and nonlinear terms. One of the most 
successful versions is arguably the "Boltzmann-Ginzburg-Landau" (BGL) framework 
recently put forward by some of us [211 l2H], where, in the spirit of weakly nonlinear 
analysis, one performs well-controlled expansions in the vicinity of ordering transitions. 
Kinetic approaches alone thus yield good deterministic "mean-field" equations but one 
still need to "reintroduce" fiuctuations in order to get bona fide mesoscopic descriptions. 

In this work, we show how this complete program can be achieved for the case of 
active nematics, i.e. systems where particles are energized individually but not really 
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self-propelled, moving along the axis of the nematic degree of freedom they carry, with 
equal probability forward or back. (Think of shaken apolar rods aligning by inelastic 
collisions \L2\-) Starting from the Vicsek-style model for active nematics introduced in 
|26j . we formulate a version of the BGL scheme mentioned above adapted to problems 
dominated by diffusion, derive the corresponding hydrodynamic equations, and study 
their homogeneous and inhomogeneous solutions. In a last section, we show how these 
equations can be complemented by appropriate noise terms using a direct coarse-graining 
approach. 



2. Kinetic approach 



2.1. Microscopic dynamics 

We consider the microscopic model for active nematics of in two space dimensions. 
This Vicsek-style model can be thought of as a minimal model for a single layer of 
vibrated granular rods [12] although it does not deal explicitly with any volume exclusion 
forces. Here, rather, pointwise particles j = 1, . . . , are characterized by their position 
X* and an axial direction 6j G [— 7r/2, 7r/2]. They interact synchronously with all 
neighboring particles situated within distance ro in a characteristic driven-overdamped 
dynamics implemented at discrete timesteps At: 



e 



t+At 



1 



-Arg 



E 



,t+At 



x* + do k] n* 



(1) 
(2) 



where Vj is the neighborhood of particle j, d^ < ro is the elementary displacement, 
n*- = (cos sin 6'*) is the nematic director, and ip and k are two white noises: 
the random angle familiar of Vicsek-style models, is drawn from a symmetric 
distribution Prfi'ip) of variance rj'^, and the zero average bimodal noise = ±1 
determines the actual orientation of motion. Both noises are delta correlated, namely 

(44') ~ (^/-jvD ~ Stt'Sjk- 

In the following, we adopt the convention [nfi]^^ = fiafif^ and label coordinates by 
greek indices, a, (3, . . . = 1,2, summing over repeated indices. 



2.2. Timescales and lengthscales 

We consider low density systems in which particles, at a given time, are either non- 
interacting or involved in a binary interaction. In this dilute limit we can neglect 
interactions between more than two particles. We also treat interactions as collision-like 
events, with the mean intercollision time 



Mesoscopic theory for fluctuating active nematics 



4 



where po is the global particle density and is shortest microscopic timescale of the 
dynamics, associated to the inversion of the rods direction of motion ~ At. This 
inter-collision time is much larger than the collision timescale 

2 

(4) 

For driven granular rods, Td may be thought of as the inverse of the shaking frequency, 
and for typical parameters it is much smaller than both the collision (tcou) and the mean 
intercoUision (rfree) timescales; at low enough densities ^ Tcow ^ Tfree- Note that the 
timescales (|4])-([3]) are different from the ones characteristic of ballistic dynamics [20] . 

To develop a kinetic approach we consider a mesoscopic timescale tb such that 
TcoU ^ ^ Tfrco- As a cousequeucc, we will treat the inversion of the direction 
of motion as a noisy term through Ito stochastic calculus [27]. We also consider a 
mesoscopic coarse-graining lengthscale ^B which, while being much smaller than the 
system size L, is larger than the microscopic scales, such as the step-size do, the mean 
interparticle distance Po and the interaction range tq. To summarize, in a dilute 
system one has 

..«r.(|)%<r««^ (5) 

and 

rfo < '^o < -^is <.L (6) 

where L is the system size and we have made explicit the condition that the typical 
coarse-graining lengthscale is is such that many particles are contained in a box of 
linear size Eb, that is poi^ ^ 1. 




2.3. Master equation 

We now write down a Boltzmann-like master equation in terms of the single particle 
probability distribution f{x,6,t), with — | < ^ < |, evolving over the timescale At ^ 
Tb- The minimal spatial resolution is such that many particles are contained in a spatial 
volume (Px centered around the position x. Moreover, we consider a dilute system, so 
that interactions (collisions) between particles are sufficiently rare to justify i) binary 
interactions (as explained above, particles then either self-diffuse or experience noisy 
binary, collision-like interactions), ii) decorrelation of the orientation between successive 
binary coUisions of the same pair of particles, that is /2(x, 9i, 02, t) ^ /(x, 6'i, t)/(x, 62, t). 
We first omit collisions and angular diffusion, only considering Eq. ([2]) to get 

/(x, e,t + /\t) = ]^ [/(x + n(e)rfo, e, t) + fix - h{9)do, 9, t)] , (7) 

where we have considered that a particle moves along one of the two orientations of n 
with equal probability. On the mesoscopic timescale ^ ~ At, Ito calculus |2Z] to 
second order gives 

dtf{x,e,t) = Dod^d^[n^{e)Me)f{x,e,t)] (8) 
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where 



(9) 



is the microscopic diffusion parameter. 



To account for angular diffusion and binary colhsions, the appropriate integrals 
need to be added to the right hand side of Eq. (|8| , 

dtf{^,e,t) = Dod^dp[h^{9)hp{e)f{x,e,t)] + h^[f] + ho,,[fj]. (lo) 

The diffusion integral describes self-diffusion which takes place at a rate \ = 1/rd 

/7r/2 POO 
do'fie') dCPi05AO'-o + (11) 
-n/2 J-oo 

where we used the simplified notation f{6) = f{x,6,t), 5-,^ is a generalized Dirac delta 
imposing that the argument is equal to zero modulo vr and -P(C) is a symmetric noise 
distribution of variance o"^, corresponding to the effective noise arising at the timescale 
tb from the sum of the microscopic stochastic contributions to angular dynamics. 
Binary collisions are described by 

/coii[/, /] = - f{e) r^'de'f{e')K{e, e') (12) 

J-w/2 



7r/2 



/2 
tt/2 



'1t/2 J-Tr/2 J ~oo 

where, for the sake of simplicity, we have used the same noise distribution P(C) as in the 
self-diffusion integral, and the out-coming angle ^ from deterministic binary collisions 
is, for -f < 61,62 < f, 

]^{6i + 62) + h{6^-62) with h{6) = i ° if 1^1 < 



,60 



if - < 



2 

< TT 



(13) 



Note that the role of the function h{6) is to ensure that ^(^1,^2) is 7r-periodic with 
respect to 61 and 62 independently. The collision kernel K{6i,62), i.e. the number of 
collisions per unit time and volume, is calculated as follows. Consider two particles 
with nematic axis n(^) and n(^') located in the volume rf^x centered around position 
X. In the reference frame of the first particle the second one diffuses either along the 
|n(^) — n(^')l oi' l^i(^) + ^{^')\ nematic axis. In unit time, taking into account the 
characteristic timescales and step-size do of its motion, it sweeps a surface (its cross 
section, which is conserved going back to the lab reference frame) equal to 

roc/o 



K{6,6') 



2ao 



\hi6) 
6 

sin — 



n 



6' 



HO) 

6-6' 



n 



cos 



where we have introduced the microscopic collision parameter 
roc/o 



"0 



(14) 



(15) 



Mesoscopic theory for fluctuating active nematics 



6 



Note that K{e, 9') = K{0 - 9') is an even function of the difference [9 - 9'), and fulfills 
the nematic symmetry, being invariant under rotation of either angle by tt. 

Before proceeding to derive hydrodynamic equations, we simplify all notations by 
rescaling time i = \t = t/rd and space x = ^x. As in [25] we also set the collision 
surface S = ^tqcIq to 1 by a global rescaling of the one-particle probability density /, 
without loss of generality. This amounts to set Aq = 1, -Do = 1 and 2q;o = 1, so that, 
dropping the tildes, our Boltzmann-like master equation now depends only on the global 
density po and the noise intensity a. 



2.4- Hydrodynamic description 



In two spatial dimensions, hydrodynamic fields can be obtained by expanding the single 
particle probability density / in Fourier series of its angular variable 9 E [— 7r/2, 7r/2|§] 

^ k=oo 

f{x,9,t) = - fkMe-'''' (16) 

TT ^ — ^ 



k=—oc 



and 



/fe(x,t) 



7r/2 



Tv/2 



d9f{x,9,t)e 



i2ke 



(17) 



The number density and the density-weighted nematic tensor field w = pQ are then 
given by 

-77/2 



p(x,t) 



7r/2 



d9f{x,9,t) = fo{x,t) 



and 



Wn{x,t) = -W22(x,t) 



Wl2(x,t) = W2l(x,t) 



d9f{x,9,t) cos(2e) = jRe/i(x,t) 
^ ^ ' d9f{x,9,t)sm{29) = hmM^,t) 

-7r/2 ^ 



(18) 

(19) 
(20) 



Note that when Im/i = the nematic field is aligned either along the x (Re/i > 0) or 
the y (Re/i < 0) axis. 

Injecting the Fourier expansion (16) in the master equation (10), one gets, after 
some lengthy calculations detailed in [Appendix A the infinite hierarchy: 



djk{x,t) = ^A/fc(x,t) + ^ ( V*Vfc+i + VVfc-i ) + 



Pfc-1 



+ -y'/g(x,t)/fc-g(x,t) 
71 ^ — ' 



PkJl 



k'Jk,q 



16g2 



(21) 



§ These A:-modes are equivalent to even harmonics if one would define particles orientation in [— tt, tt] 
in spite of the symmetry under rotations by tt (with odd ones being zero by symmetry). 
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where Pk is the Fourier transform of the noise distribution P(C) (namely, Pk 
dCP{C)e''^^) and 

J _ l + 2V2(2g-A;)(-l)gsin(f) 

''"^ ~ 1 - 4(2g - ^ 

and we have introduced the following "complex" operators 

V =dx + idy 
V* =d^- idy 
A = VV* 



22) 



The equation at order k 
1 



is thus expressed in the simple form 



dtp 



2Ap+-Re 



(23) 



and is nothing but the continuity equation for diffusive active matter with local 
anisotropy characterized by /i. 

Eq. (21 ) possesses a trivial, isotropic and homogeneous solution: p(x, t) = /o(x, t) = 
Po and fk{x,t) = for \k\ > 0. We are interested in a nematically ordered homogeneous 
solution which could eventually arise following some instability of the isotropic solution 
above. In analogy to the scaling ansatz used for polar particles [201 ES] , the interaction 



term in Eq. (21) suggests a simple scaling ansatz to close the infinite hierarchy of 



equations on fk{x,t): Near an instability threshold with continuous onset, Fourier 
coefficients should scale as fk{x,t) ~ e'*^' where e is a small parameter characterizing 



the distance to threshold. Moreover, the curvature induced current (last term of (23)) 
also induces an order e variation in the density field, p(x, t) — po ~ e- Then, assuming 



spatial derivatives to be of order e, the request that all terms in Eq. (23) are of the 



same order also fixes the diffusive structure of the scaling of time and spatial gradients: 



dt 



A 



~ e 



Using the above scaling ansatz, we proceed by discarding all terms appearing in 
(21) of order higher than e^. For /c = 1, 2 we get: 

1, 



and 



dj, = ^aA + ^VV + ai(p)/i + 6i A*/2 



= -VVi-a2(p)/2 + &2/i/i 



where the coefficients are 



ai(p) 



37r 



(2^2- 1) A - I 



p-(i-A), 



3157r 



13-9A(1 + 6V2) 



(24) 
(25) 

(26) 
(27) 
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TT 



1 

15 



(2J 



(29) 



Eq. (25) shows that at this order /2 is enslaved to /i (given that 02 > 0) and, further 

1 



a2(Po)/2 



rV^/i + 62/1/1 



(30) 



where the coefficient 02 is evaluated at the mean density po, since the 6p = p — po 



corrections are of higher order. By substituting Eq. (25) into (24) one finally gets, 
neglecting the term /^ V^/i ~ e"^, 

2' 



dtf: 



/: 



where we have introduced the transport coefficients 

71 



P 
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3^ 
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I Pi 



p 



I -Pi 



+ P2 



(1 + 6^2) A 



13 

~9 



with u 




Po + 1 - 



(31) 

(32) 
(33) 

(34) 



Note that the coefficient ^ is only a function of the average density po, as space and 
time dependent corrections are of order e^. Note also that the coefficients p and ^ are 
exactly the same as those found for the nematic field equation of nematically-aligning 
polar particles 

Eqs. (23) and (31) can be expressed in tensorial notation. To this aim, we 
introduce the linear differential operator F, such that Fn = — F22 = ^i^i — 92^2 and 
F12 = F21 = 28182, and the Frobenius inner product A : B = Aaj3Ba/3 (note that 



w : w 



w 



P and F : w = 28a8i3Wap)- After some manipulation of the terms and the 
use of Eqs. ( TOpO ), we obtain the hydrodynamic equations for the density and nematic 
field 



8tp 
8tw 



1. 1/^ 
-Ap+-(F:w) 



1 



1. 



pw — 2,^ w (w : w) H — Aw H — Fp 



(35) 
(36) 



Although the tensorial notation might be more familiar to some readers, it is in fact 
easier here to continue manipulating the complex field /i and the complex operators 

II Note that in |25j . the equations obtained are not entirely correct: (i) there is a sign error and a 
misplaced factor tt in the expression of ^; (ii) the term |V^/2 should read |A/2, where A is the 
Laplacian. In addition, let us emphasize that the Fourier coefficients Pk have a different definition in 
pS] . due to the absence of global nematic symmetry: P^. here corresponds to P2k in [53], leading to 
(only apparent) differences. 
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defined above. Moreover, in the following we drop the " " " superscript to ease notations. 
Eqs. (35) and (36) are also derived from an apolar Vicsek-style model in |29j . 



The parameter- free character of the Laplacian term in (36) means, consistent with 
our expansion in e, that the nematic phase of our system will be characterized by a 
single Frank constant [2H]. The nonlinearities studied in ^0] are therefore also absent 
to this order. The last term in Eq. (23) (or Eq. (35)), i.e. ^Re(V^/i) (or \{T : w)). 



is a curvature induced current which couples the density and the nematic field. While 
its existence was first deduced from general principles [E], here we have computed it 
directly from microscopic dynamics. Our calculations also give an exact expression for 



the corresponding transport coefficient, which is equal to the diffusive one (in Eq. (23) 



or Eq. (35)), here set to 1/2 by our rescaling. In Appendix B, we show explicitly that 



this curvature-induced current originates from the coupling of orientation with motility. 

We note finally that Eqs. ( 35p6 ) are similar to those found by Baskaran and 
Marchetti [22] but simpler, largely due to our simpler starting point. 



2.5. Homogeneous solutions 

From now on, we use for P(C) a centered Gaussian distribution of variance cx^, in which 
case Pfc = e~'^^^"^ . The linear stability with respect to homogeneous perturbations of 
the disordered solution p{x,t) = po, /i(x, t) = is given by the sign of /i(po) which 
yields the basic transition line 



\ 



In 



8{2V2- l)po + 37r 
56po + 157r 



(37) 



Note that in the dilute limit po ^ 1, where the equations have been derived, one has 
c^t ~ a/Po- 

For cr < (Jt, p > 0, and the homogeneous nematically ordered solution 

f (38) 

exists and is stable w.r.t. homogeneous perturbations. The critical line is shown in 
Fig. [T^ (black solid line). Note that for a < at, all transport coefficients (32 34) are 
positive. This will be useful in the rest of the paper. 



l/i 



3. Linear stability analysis 

We now study the linear stability of the above homogeneous solutions w.r.t. to arbitrary 



perturbations. Linearizing Eqs. (23) and (31) around a homogeneous solution, fi = 
/i,o + Sfi and p = po + Sp, one has 

dt6p = J A(5p + ^Re (V*^(5/i) (39) 



dtSfi = (po - e l/i.oH 5/i + /x7i,o Sp - 2Uifl Re (A*o + \^^5p + ^A5/i (40) 
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where /xq = /^(po) and /i' is the derivative of yU w.r.t. p. We then introduce the real and 
imaginary parts of the order parameter perturbation, 6fi = 6f[^^ + iSf[^\ and express 
the spatial dependence of all perturbation fields in Fourier space, with a wavevector 
q = {qx, qy), by introducing the ansatz 

5p{x,t) =5pqe^*+*'^^ (41) 
5/}^) (x, t) = 5/( J 6^*+^^^^ , 6fi'^ (x, t) = 6fi'l e^'^'^^ . (42) 

The stability of the stationary solution fi^ is then ruled by the real part of the growth 
rate s. 



3.1. Stability of the disordered isotropic solution 

We first study the stability of the disordered solution fi^ 
Substituting Eqs. gl]), (Q in Eqs. (Q, Q, one has 

^2 



0, in the case /iq < 0. 



q ' 



(43) 



Qy)Spci + PO 



sSf[ 



where q"^ 
simplicity qx 



(I) 
i,q 

2 



1. 2 



-qxqySpq + [ po- 



q' 



qi + qy- directions of the wavevector q being equivalent, we choose for 
q and qy = 0. From Eq. (43), one then sees that the component 



becomes independent from 5pq and 5/{ q , yielding the negative eigenvalue s = po — y- 
The eigenvalues of the remaining 2x2 block of the stability matrix are solutions of the 
second order polynomial 



+ s [g^ - po] 



4 



Po 



(44) 



In the disordered state po < 0, so that f3i and (32 are positive and one always has 
Re(s) < 0. Therefore, the homogeneous disordered solution is stable w.r.t. to all 
perturbations if po <0, i.e. a > cTf 



3.2. Stability of the ordered solution 

To study the stability of the anisotropic ordered solution, it is convenient to choose a 
reference frame in which order is along one of the axes: 

Re(/i,o) = ±y^, Im(/i,o)=0. (45) 

This solution is aligned along x, if fi^ is positive, or along y if negative. For simplicity 
we will concentrate further on the case /i o > 0, i.e., on the nematic solution aligned 
along the x axis. The real part Sf[^^ of the nematic field perturbation describes changes 
in the modulus |/i,o|, and the imaginary part 6f[^^ describes perturbations perpendicular 
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to the nematic orientation. The ansatz (41), (42) then yields the three coupled linear 
equations 

.2 1 



i,q ' 



(46) 



p7: 



1,0 - ^(Q'x ^ % 



^Pq - 



2/io + 



5/} 



(R) 

.q 



2 " 2 

We performed a full numerical stability analysis of these equations. The results are 
presented in Fig. [TJ The transition to the homogeneous solution is given by the line at- 
This solution is unstable to finite wavelength transversal perturbations of angle 16*1 > f 
between the lines at and cr^ (dotted purple line in Fig. [T|, but is stable deeper in the 
ordered phase. 

Two remarks are in order. First, the angle of the most unstable mode is here 
always perfectly |. It is thus possible to obtain the "restabilization" line ag analytically 
as shown below. Second, there is no spurious instability at low noise and/or high 
density (although we have found that such an instability appears if the truncation of 
the equations is made to the fourth order). 

To obtain the analytic expression of the line as, we write the wavevector in terms 



of its modulus q and its angle 9q, so that ql — 



g^cos2^q and 2qxqy = g^sin2^q. 
We can then analyze Eqs. (46) in the longitudinal and perpendicular wavedirections 
6'q = 0,±|, where the imaginary perturbation 5/i^q decouples from the other two. 



±- 

w , _i_ 2 , 

The latter is stable towards long wavelength perturbations, since the corresponding 
eigenvalue s = —q^jl is negative. The stability towards density and real perturbations 
depends on a 2 x 2 matrix which yields the quadratic eigenvalue equation 



+ [2/io + g^] s + 
whose solutions are 



(1 +^ 







(47) 



2fio-q^±\UfilT2fi'fi,oq^ + 



(4J 



The sign ± in front of the /i'/i,o term in Eq. (47) corresponds to the case ^q = (positive 
sign) and 6'q = | (negative sign) respectively. Note that fi' is strictly positive, as typical 
for all active matter system with metric interactions, where the interaction rate grows 



with local density. Also /io is positive and of order (see Eq. (45)). It it thus easy to 
see that in the case of large q, ^[s] < 0. For small values of q, we perform an expansion 
to order of the largest growth rate s+, obtained by taking the positive sign in front 



of the square root in Eq. (48), leading to 

,2 



2 



p 

2/io 



/i, 







(49) 
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Figure 1. (a) Basic stability diagram. The line at (solid, black) marks the linear 
instability of the disordered homogeneous solution. The ordered homogeneous solution 
is linearly unstable to large wavelengths between the at and as (dotted, purple) lines, 
and linearly stable below the as line. The (Tmin and (Tmax lines mark the domain of 
existence of the band solution ( [62| . (b) Density and order profile of the band solution 
for po = 1, a — 0.265, L = 1000; note that the lower and upper levels (pgas and pband) 
are respectively lower than pt and higher than ps, i.e. such that the corresponding 
homogeneous solution are lineally stable, (c): properties of the band solutions for 
Pq = 1: left: values of pgas (long dash, dark blue line) and pband (dashed, red line) as 
a varies between (Tmin and CTmax; right: corresponding variation of the surface fraction 

UJ. 



We can then conclude that for longitudinal perturbations (^q = 0, negative sign in front 
of fi'), the homogenous solution is stable confirming the results of numerical analysis. 
In the case of transversal perturbations (6'q = ±71/2), the stability condition is given by 

/io > ^ (50) 

meaning that close to the instability threshold of the disordered solution, when is 
positive but small, the state of homogeneous order is unstable with respect to long 
wavelength perturbations. This instability was first identified in a kinetic-equation 
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analysis by Shi and Ma [23j. Note that condition (50) is vahd up to the third order in 
e (or, equivalently, in the order parameter ||w||). It yields the stability line 

where /i2 = Ai(p = 0), = (1/0' and ^2 = (V^Po = 0)). We do not provide here 
the explicit analytical expression for as because this requires solving a sixth order 
polynomial. 

We remark that the near-threshold instability discussed above is rather generic and 
appears in "dry" active matter systems with metric interactions, as opposed to systems 
with metric-free ones, where the interaction rate is density-independent, and /i' = 
[311 12H 132] . In this case (topological active nematics), stability would be enforced by 
the positive higher order corrections fioq"^ which dominates arbitrarily close to threshold. 



4. Inhomogeneous solution 



We now show how a spatially-inhomogeneous stationary "band" solution to our 
hydrodynamic equations can be found. First we remark that our equation for the 
nematic field Eq. (31) is formally the same as that derived in [25j for polar particles 



with nematic alignment when the polar field is set to zero, as it is imposed here by 
the complete nematic symmetry of our system. We thus expect an ordered band 
solution made of two fronts connecting a linearly stable homogeneous disordered state 
(p = pgas < Pt) and a linearly stable homogeneous ordered state (p = pband > Ps) (see 
Fig. [1]) . Following [2S] , we rewrite 

p(p) = p'(p - Pt) , (52) 

with Pt = (1 — A)/p', suppose that the nematic field is aligned along one of the axes 
and varies only along y. In other words: 

Re(/i) = /i(y) , lm(/i) = , p = piy). (53) 



Eqs. (35) then becomes 



dip 



(54) 



which can be integrated to give 

p = fi + Ay + Pgas (55) 

where A and Pgas are integration constants. Furthermore, to keep the fields finite for 
\y\ —7- 00, one has A = 0. By substituting Eqs. (54) and (55) into Eq. (31) one gets 

Pt)/i-4p7i +4e/f 



dyyfl = -4:fi' (p. 



gas 



(56) 



We multiply equation (56) by dyfi and integrate it once to obtain 

I {dyhf = -2p' (Pgas - Pt) f! - tp'fl + e/f. 



(57) 
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Separating the variables we obtain 



dy = ±l '^^^ (58) 



Integration of this equation under the condition \iv!iy^±oofi{y)=^ gives after 
simphfications 



h (y) = ^i^' ''^'^ (59) 



1 + a cosh (^2?/a/p' (pt - Pga 



where a = y 1 ~ ^ (Pt — Pgas)- We still need to obtain the value of pgas which is fixed 
by the condition p (y) dy = pqL, where L is the length of the box. In the integral on 
the l.h.s we can neglect the exponentially decaying tails and integrate instead on the 
infinite domain. Furthermore, in the limit L — )■ oo we can neglect the exponentially 
weak dependence of pgas on L everywhere except the a term. We then obtain 

P,^ « ft - 1^ (1 - . (SO) 

A- -¥l^fl + ^^(p«-.0V (61) 



9v^ V V 



Substituting it back into Eq. (59) we get, under the assumption L ^ oo: 



xband 9/;' 

l + 2e-'i^cosh(y^n 3e 



and we finally obtain the ordered solution density 

Pband = A'""'' + Pgas = Pt + ^ + 26"^^) (63) 

with, as expected, pband > Pt > Pgas, which guarantees the stability of both the ordered 
and disordered parts of the solution. Note that since f^^'^'^ > the nematic order is 
parallel to the x direction (i.e. along the band orientation). This is the opposite of what 
happens in the Vicsek model, where bands extend transversally with respect to their 
polarization [15] . 

We can introduce the band fraction Q which indicates the fraction of the box 
occupied by the band. If we suppose that the front width is negligible (once again 
justified in the limit L — )■ oo), this band fraction is determined by the equation 

(Pband - Pgas) + Pgas = Po (64) 

Substituting inside the values of pgas and pband, we obtain 
9e (po-pt) + 2p^ 

n = (65) 

The condition < i7 < 1 gives us the lower and upper CTmax hmits of the existence 
of bands. As found for polar particles aligning nematically, these limits of existence 
of the band solution extend beyond the region of linear instability of the homogeneous 
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ordered solution (given by a G [dsjCrt], see Fig. [T]). In Fig. [T| we provide a graphical 
illustration of the shape and properties of the band solution. 

An important problem left for future work is the linear stability analysis of the band 
solution in two space dimensions. This is all the more important as the unpublished 
work of Shi and Ma [22] suggests the existence of some instability mechanism. 



5. Langevin formulation 

Being based on a master equation, the derivation we have discussed in the previous 
sections leads to a set of deterministic PDEs. This is a standard approach in equilibrium 
statistical physics, where the microscopic fluctuations are integrated out in the coarse 
graining process implicit in the definition of a mesoscopic cell size Ib- Fluctuations, 
when needed, can be eventually introduced as an additive, delta correlated stochastic 
term as in Ref . [16] . However, the presence of large density fluctuations [IT] suggests 
that fluctuations may not be faithfully accounted for by some additive noise term. The 
precise nature of noise correlations at the mesoscopic level cannot be safely overlooked 
in non-equilibrium systems, as it is known that stochastic terms multiplicative in the 
relevant fields can radically alter the universality class of mesoscopic theories [53] . 

In this section, we perform a direct coarse-graining of the microscopic dynamics in 
order to compute the (multiplicative) stochastic terms which emerge at the mesoscopic 
level. We however restrict the computation to the stochastic terms emerging from 
the coUisionless dynamics. For real-space coarse-graining, we make use of a smooth, 
isotropic, normalized (to one) filter (7s (r) decaying exponentially or faster for r > s, e.g., 
a Gaussian of width s. The fluctuating coarse-grained density and nematic order field 
are then defined as 

N 

p{x,t)^J29M-^) (66) 



and 

N 

w(x, ■) 

i=l 

where we have introduced the microscopic traceless tensor 



,t)^ J](7,(x*-x)Q*, (67) 



- "^"^ "2-2 sin 26^ - cos 26^ ' = ' ^^^^ 



5.1. Density field fluctuations 



The correlations of density field fluctuations can be derived by generalizing an approach 
first outlined by Dean [3l] for Brownian particles. As mentioned above, we use the 
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collisionless dynamics. We are interested in the time evolution of the density field (66), 
which is given by 

N N 

p(x, t + At) = ^ ^,(x*+^* - ^) = E 9si< + Ax* - x) (69) 



1=1 



i=l 



,t+At / 



where Ax* 

Expanding up to second order in powers of Ax* according to Ito calculus [2Zj and by 
virtue of Eq. ^ one has 

dtp{x,t)=To{x,t)+Ti{x,t), (70) 

where 



"o 



N 



= ^ J2^hUn^\pd^dpgM - x) 

1=1 



and 



N 



T,{x,t) = ^Y.< (n*-V)^7.(x*-x) 

1=1 



(71) 



(72) 



Note that derivatives are taken w.r.t. the argument of the function Qs, and not w.r.t. x. 
The second order term Tq yields the deterministic part of the density dynamics. By 
Eqs. ( |66||67t and the definition of the microscopic nematic tensor Q [Eq. (p8j)] one easily 
gets 



To = §^(r:w) 



(73) 



that is, the right hand side of the diffusion Eq. (35). The first-order term Ti gives rise 
to the (zero average) stochastic term we are interested in. At this stage, Ti is not a 
simple function of the mesoscopic fields; however, following Ref. [Mj it is possible to 
show that its two point correlation can be recast as a function of p and w. Averaging 
over the random numbers k*, we have, in the limit s — )■ 0, 



(Ti(x,t)Ti(y,t')) = df-^^-^ ■ ■ 9s{< - x)^7.(x* - y) 



N 



dl^\^ E ■ ■ ^y) (^?.(x-y)^7.(x*-x)) . (74) 

i=l 



i=l 

N 



Using Eq. (68), one then finds, approximating the filter by a Dirac delta in the limit 
s ^ 0, 



(Ti(x,t)Ti(y,t')) = df-^^^d^dp 



(5(x-y) ( Waf}{x,t) + ]^p{x,t)5af3 



(75) 



We can rewrite the noise term Ti in the stochastically equivalent (i.e., with the same 
correlations on the mesoscopic scale) form 

Ti(x,t) = V-h(x,t) (76) 
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(77) 



where h is a Gaussian, zero-average vectorial noise, delta-correlated in time with 
correlations 

(/Z JX, t)hp{)l, t')) ^y^Kt- t') - y) (^a/3(x, t) + ^p(x, t)^ . 

Such a noise term can finally be expressed in the more convenient form 

/ia(x, t) = if„/3(x, t)h,3{x, t) , (78) 

where the Gaussian noise h has correlations independent from the hydrodynamic fields 

Ch^{x,t)hp{x',t')) = 2DoSo.p6{t-t')6{x-x') (79) 

and the tensor K is implicitly defined from the relation K-K = (p/2)I + w (with I 
being the identity matrix). In the limit of small w considered here, we can expand K 
to first order in w, yielding 

K = ^p'/'(l+-). (80) 



72^ V P 

The divergence term V- appearing in Ti refiects global density conservation, while the 
proportionality of noise variance to number density can be interpreted as a consequence 
of the central limit theorem. Adding up the two contributions, one finally gets 

9,p=^(r:w) + ^VV + V- (K-h). (81) 
5.2. Nematic field fluctuations 



We next discuss fiuctuations of the nematic tensor. As seen from Eq. (67), w is a 
function of the 2N microscopic stochastic variables x- and — through the microscopic 
nematic tensor (68) — 6**, whose dynamics is given by Eqs. ([T])-(|2|. According to Ito 
calculus, one has 

dtw = rio + rii + ris (82) 

where flo is the deterministic part of the coarse-grained collisionless dynamics (which 
we do not write here explicitly), while fii and 0,2 are two stochastic contributions, 

2 ^ 

where ipi and are the microscopic noises and 

A=| ° 1 • (85) 
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Note that in fii we have retained only the linear contribution in the microscopic noise 
ipj. We first focus on the stochastic terms fii. On coarse-graining scales, averaging over 
the microscopic noise ipj, correlations of fii are given by 

r)(f-f') ^ r 1 r 



Td 



i=l 



7(5 



N 



V^^^.(y-x)5^^.(x*-x) 



Td 



i=l 



A-Q* A-Q* (86) 



7<5 



To evaluate this correlator, we determine the average value {J2i9s (A ■ Qf){-^ ■ Q*)), in 
the framework of the deterministic dynamics studied in Sect. [2| namely 



N 



^c/.(x--x) 



i=l 



A Q* 



J al3 



A Q* 



7(5 



d9f{x,9,t) [A-Q(^)^ [A-Q(^^)]^,(87) 



After some rather lengthy calculations, using the closure equations (30), one finds 

2b 

d9f{x,9,t) [A ■ Q(6')]„o [A • Q{9)] . = pJafSyS H ^ [iw^^Wf,^)Ja(3^5 - 2Wal3W^5] 

02 



+ -. [^fiuW^uJalS-yS — ^a/SW^S — ^-ySWa/s] ,(88) 

4a2 



where we have introduced the tensor 



(89) 



which plays the role of a unit tensor for the double contraction of symmetric traceless 
tensors, e.g., Wai3 = JaiSfiuWf^,^. In order to characterize the noise fii, we introduce the 
following change of variables: 

[rii(x, t)]^^ = H^p^,{x, t) Vt^,{x, t) (90) 
where f2 is a tensor ial symmetric traceless white noise, such that 

(fi„^(x, t)fi-,5(y, t')) = 2D5{x - y) 5{t - t') J^^^s, (91) 
with D = 27]"^ /Td- The correlation of fii then reads 

([f^i(x, t)]„/3 [f^i(y, t%s) = 2D6{x - y) 6{t - t') H^p^,{x, t)H,s,A^, t) (92) 



By identification with Eq. (86), and using Eq. (88), one eventually finds for H 



0.2 P 



1/2 



5^2 P 



Q/3J 



(93) 



Note that, in agreement with the central limit theorem, fii is (at least to first order in 
w) proportional to the square root of local density. 

The second stochastic term f22, finally, can be treated similarly, but it would give 
rise to a conserved noise (due to the presence of V terms) akin to the one discussed for 
the density equations, thus related to density fluctuations affecting the w = pQ field. 
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We discard such conserved term as irrelevant (in the renormahzation group sense) with 
respect to the non-conserved multiphcative noise fii. 

In order to write down the complete Langevin equation, one also needs to evaluate 
the contribution of the deterministic part fio- However, expressing this contribution in 
terms of the fluctuating fields p and w turns out to be a very complicated task. One 
should also take into account collisions between particles, and not only the coUisionless 
dynamics described by fio- Then some further approximations would be required to 
treat the non-linear part of the dynamics. 

In addition, microscopic collisions could provide a further fluctuation source due 
to disorder below the coarse-graining scale. While we conjecture them to be irrelevant, 
we leave a final settlement of this difficult problem for future work, and use for the 



deterministic part of the dynamics the terms the hydrodynamic equation (36), derived 
from the Boltzmann approach. 

We thus finally obtain the stochastic equation for the nematic field 

dt-w = /iw - 2^ w (w : w) + -Aw + -Tp + H : i7 . (94) 

2 8 

A few remarks are in order: first, our expressions of the noise amplitudes K and H 
(Eqs.([80)) and ([93])) suggest that the stochastic terms might be better expressed in 



terms of the field Q, rather than w = pQ; second, Eqs. (81) and (94) are also derived 
from an apolar Vicsek-style model in |29j . 

In spite of the limitations listed above, the present approach already provides us 
with useful information on the statistics of the noise terms, which is seen to differ 
significantly from the white noise postulated on a phenomenological basis in previous 
works. On top of the overal p^/^ dependency, our calculation reveals a non-trivial 



dependence of the correlation of the noise on the nematic order parameter [see Eqs. (80 



81, 92, 93) 



6. Conclusions 

To summarize, using as a starting point the simple active nematics model of ^26j, we 
have demonstrated how one can derive in a systematic manner a continuous mesoscopic 
description: We formulated a version of the Boltzmann-Ginzburg-Landau approach put 
forward in [211 [25] for this case where (anisotropic) diffusion dominates, deriving a simple 



hydrodynamic equation for the nematic ordering field -Eq. (36). We have then used 
a direct coarse-graining approach to endow the hydrodynamic equations with proper 
noise terms. 

The next stage, left for future work, consists in studying the stochastic PDEs 
obtained. At the linear level, it is clear that in the long wavelength limit, standard 
results on giant density fluctuations ^/Tj are recovered. However, the large amplitude of 
density fluctuations calls for a non- linear analysis (which turns out to be very difficult), 
where the density dependence of the noise derived in Sect. [5] may play an important 
role. Ideally, one should try to tackle this issue by applying methods from field theory 
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and renormalization group analysis. In addition, we note that the multiphcative nature 
of the noise may also affect finite-wavelength properties, like coarsening behavior. The 
analysis of the stochastic PDEs can be done numerically, but some care must be taken 



when dealing with the multiplicative, conserved noise terms in (81). 

Pending such attempts, some remarks and comments are already in order: like all 
previous cases studied before, the hydrodynamic equations found exhibit a domain of 
linear instability of the homogeneous ordered solution bordering the basic transition line 
cTf This solution does become linearly stable deeper in the ordered phase (for a below 
(Js). Moreover, we have found that the long wavelength instability of the homogeneous 



ordered solution leads to a nonlinear, inhomogeneous band solution -see Eq. (62)- and 
that this band solution exists beyond the [us, at] interval. These coexistence regions 
suggest, at the fluctuating level, discontinuous transitions. 

This seems to be at odds with the reported behavior of the original microscopic 
model: (i) the order/disorder transition has been reported to be of the Kosterlitz- 
Thouless type [26]; (ii) there is no trace, at the microscopic level, of the existence of a 
non-segregated, homogeneous phase; (iii) coming back to giant number fluctuations, we 
note that the standard calculation is made in the homogeneous ordered phase whereas 
the numerical evidence for them reported in [26] appears now to have been obtained in 
the inhomogeneous phase. All this calls for revisiting the simple particle-based model 
and, eventually, understanding its behaviour in the context of the stochastic continuum 
theory constructed here. 
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Appendix A. Fourier expansion of the master equation 

We provide in this Appendix details of the Fourier expansion of the master equation 
(10), leading to Eq. (21). Multiplying Eq. (10) by e*^^ and integrating over 6, one gets 



tt/2 



dtfk = d^dp I dee''''K{e)np{e)f{x,e,t) 

-Tr/2 



n/2 



tt/2 



Tr/2 



+ / de e''''hm[f] + / de e^'^'lconif, f] ■ {A.l 



Tr/2 



In the following, we successively compute each term of the r.h.s. of Eq. (A.l). 
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Appendix A. 1. Diffusion-like term 
Let us define Qapi^) as 

We then have 



Qn{0,t) = -Q22{e,t) = -cos2e 



Qu{0,t)=Q2i{9,t) 



sin 2^ 



g«2e _ g-j29 



4i 



As a result, 



tt/2 



Tr/2 



22 



(A.2) 



(A.3) 



6. 



-7t/2 J-tt/2 \ ^ 



(A.4) 



Appendix A.2. Self-diffusion term 
We have rather straightforwardly 

dOe^'^'hislf]-- 

-tt/2 



where 

/oo 
■oo 

is the Fourier transform of -P(C)- 



/7r/2 /"OO 
-n/2 J -oo 



Pk-l 



fk 



(A.5) 
(A.6) 



Appendix A.3. Binary collisions term 

Let us spht the Fourier transformed colhsion integral into an outgoing (negative) collision 
term ijr^ and an ingoing (positive) collision term Z^^''. A direct integration of the 
outgoing collision term yields, using K{9, 9') = K{6 — 9'), 

d9e''''fi9) d9'f{9')k{9-9') = -1 ^ i^JJ,_^(A.7) 



r(-) 



-tt/2 



-tt/2 



TT 



where Kg is the Fourier coefficient of K{9 — 9') given by, using Eq. (14), 

^^/2 _ „ r e-9' 9-9' 



d9e'^''^ 



tt/2 



sm 



cos 



4 



1 - 16g2 



(A.8) 



Then, the calculation of the ingoing collision term requires a few steps. After integration 
of the (generalized) Dirac delta 5^, we have 



d9, d92e'''''^'^''^^f{9,)K{9,-92)f{92) . 

-tt/2 J-tt/2 



(A.9) 
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By the change of variables (p = 0i — 02, one gets 

'•7r/2 1-^/2-92 



r(+) 



^2me2+<i>fi2)f(^Q^ + <p)K{(P)f{e2) . (A. 10) 



Pk / de2 

J-tt/2 J-7r/2-e2 

Using the vr-periodicity of the integrand with respect to 0, we can change the integration 
interval on 0, yielding 

"1^/2 i-Tv/2 



r(+) 



Pk I de2 

-■k/2 J-7r/2 



^i2fc* (6*2 +0,6*2) 



/(^2 + 0)i^(0)/(^2) 



On this interval of (j), one has from Eq. (13) 

^(^2 + 0,^2) =^^2+^ . 



Expanding / in Fourier series [see Eqs. ( 16p7)], we get 



r(+) 



Pk 
7r2 



/■k/2 /■7r/2 



-7r/2 



tt/2 



The integral over 62 is equal to v:5k,q+qi- Defining 

"71/2 



J k,g 

we finally obtain 



Jik-2q 



-7r/2 



7r 



■Jk,qfqfl 



k~q 



The coefficient J^^g can be computed explicitly, leading to 



Jk,q = 4 

Note finally that jo,g -- 



1 + 2^2(2^ - A;)(-l)« sin (^) 



1 -4(2g- A;)2 



Kq. 



(A.ll) 



(A. 12) 



J(k-2qHj^^^^ . (A. 13) 



(A.14) 



(A.15) 



(A.16) 



Appendix B. Curvature-induced current and equilibrium limit 

In this Appendix, we show explicitly that the curvature-induced current, that is the term 
|Re ( V*^/i ) appearing in the continuity equation (23), originates from the coupling of 



orientation with motility. To this aim, we consider a slightly generalized microscopic 
process w.r.t. Eqs. ([l| |2|, where particles are also allowed to move perpendicular w.r.t 
to the nematic tensor. Replace Eq. ^ by 

x*+^* = x* + rfoR(^D (B.l) 

where R(^) is a stochastic operator defining the coupling between orientation and 
particle motion, 

n(9) w.p. p/2 
—n{9) w.p. p/2 
h^{9) w.p. (l-j5)/2 
-n^(^) w.p. (l-j9)/2 



(B.2) 



Mesoscopic theory for fluctuating active nematics 



24 



where < p < 1, w.p. stands for "with probabihty" and n-'-(^) = + 7r/2) is the 
perpendicular director. The standard active nematic case is recovered for p = 1, while 
p = 1/2 corresponds to an isotropic random walk, a case for which motion is decorrelated 
from order. The corresponding collisionless master equation reads 

/(x, 0,t + At)^^ [fix - n{9)do, d, t) + /(x + h{d)d^, d, t)] 



+ 



(1-p) 



[f{x-ri'-{e)do,e,t) + f{x + h^{e)do,9,t)\ . (b.3) 



By making use of Ito calculus, one gets at the mesoscopic timescale tb 



dtf{x,9,t) = {2p-l)do.d0 



2 



/(x,^,t) + ^A/(x,^,t) 



(B.4) 



where we have used the identity n'^{0)n^{0) = Sa/3 — ha{0)hi3{6). By considering the 
zeroth-order Fourier term of / (for which collision and angular diffusion terms vanish), 
one obtains the continuity equation 



9,p=^Ap+^Re(V*Vi) 
which shows that the non-equilibrium current vanishes for p 



(B.5) 



